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The  solid  particle  solar  receiver  (SPSR)  is  a  direct  absorption  central  receiver  that  uses  solid  particles 
enclosed  in  a  cavity  to  absorb  concentrated  solar  radiation.  The  SPSR  is  a  candidate  for  applications  of 
solar  energy  in  a  thermo-chemical  water-splitting  process  to  produce  hydrogen.  This  paper  presents  a 
review  of  the  study  on  SPSRs,  including  the  idea  originality,  design  concepts,  advantages  and 
disadvantages,  the  solid  particle  identification,  a  conceptual  design  in  Sandia  National  Laboratories  and 
detailed  studies  performed  on  this  design.  The  geometry,  particle  size,  calculating  domain  selection,  the 
wind  effect,  the  aerowindow  and  other  factors  which  influence  the  cavity  efficiency  have  been  studied 
and  the  results  are  presented. 
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1.  Introduction 
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Nomenclature 

CD  drag  coefficient 

dp  particle  diameter 

J(r,s')  total  intensity 

lbx  black  body  intensity 

n  refractive  index 

Nu  Nusselt  number 

Pr  Prandtl  number  of  gas  flow 

r  position  vector 

Rep  relative  particle  Reynolds  number 

s  direction  vector 

uj  velocity  of  air  flow  (m/s) 

Uj  mean  velocity  of  air  flow  (m/s) 

u'j  fluctuant  velocity  of  air  flow  (m/s) 

uP'i  particle  velocity 

T  temperature 

f  mean  temperature 

X,Y,Z  coordinates 

Greek  symbols 

a  thermal  conductivity 

ax  spectral  absorption  coefficient 

£p  particle  emissivity 

A  wavelength 

pt  viscosity 

ptt  turbulent  viscosity 

v  kinematic  viscosity 

p  density 

pp  density  of  particle 

a  Stefan-Boltzmann  constant 

gs  scattering  coefficient 


energy  systems.  The  idea  originality  of  an  SPSR  is  to  design  an 
efficient  receiver  as  an  important  part  in  a  water-splitting  (WS) 
thermo-chemical  (TC)  process  to  produce  hydrogen,  which  is  a 
promising  energy  carrier  for  the  future  energy  supplies.  An  SPSR  is 
a  direct  absorption  central  receiver  that  uses  solid  particles 
enclosed  in  a  cavity  to  absorb  concentrated  solar  radiation.  An  SPSR 
is  an  ideal  candidate  for  applications  of  solar  energy  including  fuels 
and  chemicals  production,  Brayton  cycle  electricity  generation,  or 
industrial  process  heat  applications.  The  solid  particles  are  heated 
in  an  SPSR  and  serve  as  a  heat  transfer  and  storage  medium  or  as  a 
substrate  on  which  chemical  reaction  may  be  performed  directly  in 
solar  WSTC  processes. 

This  paper  presents  a  review  of  the  study  on  SPSRs,  including 
the  idea  originality,  the  experimental  work  and  numerical  studies, 
and  the  milestones  achievements.  Furthermore,  the  detailed 
studies  of  a  US  DOE  funded  project  on  a  conceptual  design  of  an 
SPSR  have  been  reviewed,  which  is  a  part  of  the  Solar  Thermo- 
Chemical  Hydrogen  (STCH)  project  based  on  the  cooperation  of 
Sandia  National  Laboratories  (SNL)  and  University  of  Nevada,  Las 
Vegas  (UNLV)  [  1  ].  The  review  will  include  the  conceptual  design  of 
an  SPSR  in  the  STCH  project,  the  geometry,  the  particle  properties, 
the  calculating  domain  selection,  the  wind  effect,  the  protection  of 
an  aerowindow  and  other  factors  which  influence  the  cavity 
efficiency.  The  physical  model  and  mathematical  models  are 
analyzed  and  the  solution  schemes  are  presented  with  the 
parametric  study  on  different  operating  conditions  and  model 
parameters  [2]. 


2.  Hydrogen  as  an  energy  carrier 

The  drawbacks  to  the  fossil  fuel  use  include  limited  reserves, 
and  carbon  dioxide  emissions,  which  is  a  greenhouse  gas 
responsible  for  the  global  warming  and  environmental  pollutions. 
Hydrogen  is  a  promising  and  clean  energy  carrier,  which  could 
potentially  replace  the  use  of  fossil  fuels  in  the  transportation 
sector  in  the  future.  The  US  president  George  W.  Bush’s  Hydrogen 
Fuel  Initiative  [3],  announced  in  2003,  envisions  the  transforma¬ 
tion  of  the  nation’s  transportation  fleet  from  a  near  total  reliance 
on  petroleum  based  products  to  steadily  increasing  use  of  clean¬ 
burning  hydrogen.  America’s  energy  security  can  be  dramatically 
improved  with  the  promising  hydrogen  projects  by  significantly 
reducing  the  reliance  on  imported  oil,  as  well  as  reducing 
greenhouse  gas  emissions  and  helping  to  clean  the  air. 

Life  cycles  assessment  (LCA)  of  hydrogen  production  was  first 
considered  by  Koroneos  in  Ref.  [4].  LCA  is  a  powerful  tool  to  help 
evaluate  the  impact  from  a  process  or  from  production  and  use  of  a 
product.  It  consists  of  goal  definition  and  scoping  which  defines  the 
product,  process  or  activity;  inventory  analysis  which  identifies 
material  usage  and  environmental  releases;  impact  analysis  which 
assesses  the  human  and  ecological  effects  of  energy,  water  and 
material  usage;  and  last  interpretation,  which  evaluates  the  results 
of  each  analysis.  Hydrogen  production  methods  have  been 
compared  among  natural  gas,  renewable  energy,  electrolysis, 
and  fuel  in  the  research.  The  use  of  wind,  hydropower,  and  solar 
thermal  energy  were  proved  to  be  the  most  environmentally 
friendly  and  efficient  methods  [4]. 

Hydrogen  has  much  more  advantages,  as  a  promising  future 
energy  carrier,  particularly  when  coupled  with  fuel  cells.  However, 
there  is  no  environmentally  attractive,  large-scale,  low-cost,  and 
high-efficiency  hydrogen  production  process  available  for  com¬ 
mercialization  currently  [5].  Hydrogen  can  be  manufactured  from 
a  range  of  energy  sources  such  as  fossil  fuels,  biofuels,  renewable 
energy,  nuclear  energy  via  electricity  [6-8].  Hydrogen  can  also 
readily  be  produced  from  synthesized  hydrogen  carriers  such  as 
methanol,  ammonia  and  synfuels.  The  hydrogen  production  is 
mainly  based  on  fossil  fuels  and  most  specifically  on  natural  gas  at 
present  [4].  There  are  no  environmental  benefits  if  the  hydrogen  is 
derived  from  natural  gas  or  fossil  fuel  reforming,  because  of  the 
emissions  of  C02  from  these  production  routes.  In  addition, 
hydrogen  produced  from  fossil  fuel  contains  trace  contaminants 
(CO)  that  are  poisons  for  catalysts  used  in  PEM  fuel  cells  electrodes 
[1].  A  detailed  review  of  the  renewably  driven  hydrogen  systems 
and  the  modeling  approaches  have  been  presented  by  Deshmukh 
[7].  The  use  of  wind,  hydropower  and  solar  thermal  energy  for  the 
production  of  hydrogen  are  the  most  environmental  friendly 
methods.  Among  which,  the  WSTC  cycle  is  one  of  the  most 
promising  methods  to  produce  hydrogen,  especially  when  it  is 
driven  by  solar  energy. 

3.  Solar-driven  water-splitting  thermo-chemical  (WSTC)  cycles 
for  hydrogen  producing 

Electrolysis  produces  pure  hydrogen  but  suffers  from  thermo¬ 
dynamics  inefficiencies.  The  projected  efficiency  for  the  electro¬ 
lytic  path  is  about  36%.  However,  the  current  combined  efficiency 
does  not  exceed  20-25%.  The  overall  efficiencies  of  the  most 
investigated  WSTC  cycles  are  in  the  range  of  45-50%  [1,5]. 
Therefore,  a  WSTC  cycle  is  one  of  the  most  promising  methods  to 
produce  hydrogen  from  the  point  view  of  efficiency.  In  a  WSTC 
cycle,  water  and  heat  are  the  inputs,  and  hydrogen  and  oxygen  are 
the  only  outputs.  The  other  chemical  and  reagents  are  recycled  in  a 
closed  WSTC. 

As  the  opposite  beneficial  aspect  of  the  oil  crisis,  the  interest  in 
WSTC  processes  boomed  in  the  late  1 970s  and  early  1 980s  [9-11]. 
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The  researches  in  WSTC  then  slowed  down  worldwide  till  the  last 
review  published  in  1989  [11].  Some  continuing  work  has  been 
carried  out,  but  mostly  in  those  countries  where  the  dependence 
upon  foreign  energy  sources  is  of  the  national  concern,  such  as 
Japan.  Several  hundred  WSTC  processes  have  been  proposed  for 
the  hydrogen  production.  However,  the  technical  feasibility  study 
and  efficiency  evaluation  have  been  executed  on  only  a  few  of 
these  cycles  [5]. 

Typically,  a  WSTC  process  to  produce  hydrogen  requires  a 
heat  input  in  excess  of  700  °C.  This  heat  can  be  either  supplied  by 
solar  or  by  nuclear  energy.  Historically,  processes  receiving  the 
most  attention  were  those  having  a  maximum  temperature  of 
around  900  °C  that  could  be  linked  with  a  nuclear  energy  input. 
The  TC  processes  evaluation  and  screening  were  generally 
conducted  in  the  frame  of  the  coupling  with  an  advanced 
high-temperature  nuclear  reactor  as  the  primary  energy  source 
[9-11].  Thus,  the  maximum  cycle  temperature  level  was  the 
deemed  optimum  for  an  advanced  high-temperature  nuclear 
reactor,  i.e.  about  850  °C.  Hence,  the  range  selected  for  the 
optimum  maximum  temperature  in  process  was  750-850  °C 
[12].  The  Adiabatic  UT-3  cycle  was  proposed  at  the  University  of 
Tokyo  [13,14]  and  selected  by  JAERI  for  further  development.  The 
predicted  efficiency  of  the  Adiabatic  UT-3  process  varies  between 
35%  and  50%  depending  upon  the  efficiency  of  membrane 
separators,  which  are  under  development,  and  also  depending  on 
whether  electricity  is  co-generated  along  with  the  hydrogen.  The 
UT-3  cycle  is  based  on  two  pairs  of  chemical  reactions.  The  first 
two  ensure  the  formation  of  hydrobromic  acid  releasing  oxygen, 
and  the  other  two  ensure  the  reduction  of  water  by  a  bromide. 
The  maximum  temperature  of  this  cycle  is  about  742  °C  [15], 
which  can  be  driven  by  either  the  nuclear  energy  [13,14]  or  the 
solar  energy  [16-19]. 

The  Sulfur-Iodine  (SI)  cycle  was  developed  by  General  Atomics 
[14].  The  SI  cycle  [20]  remains  the  cycle  with  the  highest  reported 
efficiency,  theoretically  52%,  based  on  an  integrated  flow-sheet 
with  suggested  process  improvements  that  could  increase  the 
efficiency  and  lower  the  capital  cost  [21].  The  SI  cycle  has  a 
maximum  reaction  temperature  of  850  °C  [22].  Therefore,  the 
primary  energy  input  has  been  designed  originally  to  be  supplied 
from  the  advanced  nuclear  power  station  (4th  generation)  [23,24]. 
The  hybrid  sulfur  (Hy-S)  processes,  proposed  by  Westinghouse 
[25],  is  another  typical  cycle  of  the  hydrogen  production  process, 
having  a  same  maximum  temperature  as  an  S-I  cycle  [26].  The  Hy-S 
cycle  is  a  variant  of  the  sulfur-iodine  cycle  and  it  involves  a 
thermochemical  and  an  electrochemical  reaction.  Actually,  the 
energy  for  SI  cycle  and  the  Hy-S  cycle  can  be  supplied  by  both  solar 
energy  and  nuclear  energy  from  their  maximum  temperature 
requirements. 

The  solar  energy  is  obviously  a  much  “cleaner”  energy 
compared  to  the  nuclear  energy  as  an  energy  input  for  the  WSTC 
process.  Moreover,  many  of  the  TC  reactions  require  a  minimum  of 
over  1000  °C  or  even  higher  than  2000  °C,  such  as  the  Zinc-Oxide 
cycle  [27].  These  kinds  of  TC  cycles  are  considered  as  solar  energy 
only  in  a  temperature  higher  than  1000  °C  is  required.  Therefore, 
WSTC  cycles  may  constitute  the  optimized  choice  for  hydrogen 
production  [28]. 

Viable  water  splitting  cycles  have  been  identified  and  evaluated 
basically  at  a  high  temperature,  which  can  be  driven  by  solar 
energy  and  can  offer  potential  for  efficient  and  large-scale 
production  of  hydrogen  [5].  Studies  about  two-step  cycles  at  a 
high  temperatures  were  related  to  CdO/Cd  (1400  °C)  [29],  ZnO/Zn 
(2000  °C)  [30-32],  Fe304/FeO  (2200  °C)  [33-36]  initially  proposed 
by  Nakamura  [37].  The  use  of  mixed  metal  oxides  (e.g.  zinc-ferrite, 
manganese-ferrite,  or  nickel-ferrite)  cycles  decreases  the  tem¬ 
perature  of  the  activation  step  while  the  reduced  oxide  remains 
reactive  for  the  H2  generation  step  [38-46].  A  project  funded  by  the 


DOE  was  initiated  in  the  USA  recently.  This  project  aims  at 
identifying  and  testing  the  high-temperature  solar-driven  TC 
cycles  [47].  A  computer  searchable  database  including  280  cycles 
[48]  was  developed  as  the  tool  to  organize  the  cycle  data.  The 
feasibility  of  the  main  part  of  the  cycles  has  not  been  proved  yet 
and  only  a  few  were  tested  experimentally.  The  basic  screening 
and  testing  work  on  the  promising  solar-driven  WSTC  cycles  can  be 
found  from  Refs  [5,49]. 

4.  Central  receivers  for  WSTC  cycles 

Solar  energy  is  much  more  environmentally  and  technically 
acceptable  as  a  heat  source  for  a  WSTC  process,  especially  at  a  high 
temperature.  However,  solar  thermal  energy  must  be  collected  by 
concentrated  methods  in  order  to  achieve  reasonable  tempera¬ 
tures  in  a  working  medium. 

The  prevailing  wisdom  with  regard  to  the  design  of  WSTC 
hydrogen  production  facilities  is  that  bigger  is  better  due  to  the 
economies  of  scale  and  operational  considerations.  Therefore,  the 
preference  of  the  WSTC  processes  will  be  those  which  are  suitable 
for  a  central  receiver  platform  [50].  A  central  receiver  facility  has 
the  potential  to  offer  thermal  storage  and  “around-the-clock” 
operation.  Concentrated  solar  energy  can  be  collected  by  a  solar 
central  receiver  system  which  uses  a  field  of  heliostats  to  reflect 
and  focus  sunlight  onto  a  receiver  located  on  top  of  a  solar  tower 
[2].  In  this  manner,  the  collected  sunlight  can  reach  hundreds  of 
suns  at  the  surface  of  the  solar  receiver.  The  concentrated  solar 
energy  is  capable  of  attaining  the  temperature  requirement  for  the 
TC  reaction  cycle  at  2000  °C  or  even  higher.  The  principal 
challenges  for  central  solar  receivers  in  a  WSTC  process  include 
the  selection  a  suitable  work  media  for  heat  transfer  and  storage 
and  the  design  and  efficiency  improvement  of  a  solar  receiver.  A 
working  medium  is  needed  to  transport  the  energy  for  the 
reactions  that  need  heat  source.  Previous  studies  on  the  working 
medium  in  a  solar  central  receiver  have  been  on  gases  or  liquids 
which  need  to  be  transported  by  tubes  or  pipes.  There  are  working 
demonstrations  of  systems  using  liquid  as  the  heat  medium  like 
the  water/steam  central  receiver  system,  located  near  Barstow, 
California,  which  produces  electrical  energy  [51].  The  other 
common  working  fluids  are  molten  salts,  molten  metals,  and  air. 
The  term  solid  particle  solar  receiver  is  another  concept  of  a  heat 
energy  transfer  medium  with  great  interest.  Sand-size  refractory 
particles  are  falling  down  freely  inside  a  solar  receiver  to  form  a 
curtain  that  directly  absorbs  the  solar  insolation.  An  SPSR  has  many 
favorable  characteristics  compared  to  other  solar  receivers.  For 
example,  the  solid  particles  absorb  the  concentrated  solar  flux 
directly  and  the  particles  themselves  serve  as  the  storage  medium 
in  an  SPSR  [2].  Therefore,  a  conceptual  simplicity  can  be  obtained 
and  a  high  temperature  (above  1000  °C)  can  be  obtained,  which  is 
not  allowed  with  other  systems,  since  a  fluid  conduit  and 
additional  piping  and  heat  exchanger  equipment  are  not  necessary 
[52].  Moreover,  there  are  other  advantages,  such  as,  the  economic- 
ity  and  abundance  of  possible  receiver  materials,  the  possibility  of 
building  a  simple  thermal  storage  system  (a  tank  filled  with  the  hot 
particles),  the  absence  of  piping  and  heat  exchange  problems  such 
as  those  related  to  the  use  of  molten  salts  (freezing  problems),  air 
or  steam  (poor  heat  transfer).  Furthermore,  there  is  no  excessive 
thermal  stress  on  the  plant  components  since  the  heat  is  directly 
absorbed  by  the  particles  without  heating  receiver  pipes  or 
cavities. 

However,  the  existing  open  aperture  lowers  the  overall 
efficiency  by  convection  heat  transfer.  It  is  necessary  to  minimize 
heat  loss  from  the  particles  to  the  air  within  the  receiver  as  this 
lowers  the  overall  efficiency  [53].  Apparently,  the  solid  particle 
receiver  addresses  both  of  the  two  principal  challenges  of  a  central 
receiver. 
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Fig.  1.  Simplified  schematics  of  the  solar  falling  particle  reactor  (FRP)  [53]. 


The  simplified  schematics  of  the  solar  falling  particle  reactor 
(FRP)  system  are  shown  in  Fig.  1  [53].  An  SPSR  is  usually  mounted 
atop  a  tower  to  capture  the  focused  solar  energy  from  a  field  of 
heliostats.  The  solid  particles  are  lifted  to  the  top  of  the  SPSR  and 
then  fall  as  a  cloud  or  curtain  through  the  receiver  where  they  are 
irradiated  directly  by  concentrated  sunlight  passing  through  the 
aperture  of  the  SPSR.  The  particles  are  typically  dark  in  color  and 
made  of  a  ceramic  material  such  as  sintered  bauxite.  The  heated 
particles  can  then  be  stored  and  run  through  a  heat  exchanger  to 
provide  thermal  energy  input  to  a  process  [2]. 

There  are  also  many  other  kinds  of  solar  receivers  under 
investigation.  In  some  designs  of  the  central  receivers,  the  heat 
absorbing  and  reaction  progress  occur  in  the  same  place 
[4,26,54,55].  For  example,  an  atmospheric-open  cyclone  solar 
receiver-reactor  has  been  designed  and  tested  by  Steinfeld  to 
decompose  the  carbonates  and  hydroxides  of  Group  IIA  elements 
in  the  periodic  table  [56].  The  receiver-reactor  is  tested  with  the 
decomposition  of  calcite  which  needs  a  reaction  temperature  of 
1170K.  The  energy  absorption  efficiency,  based  on  the  energy 
incident  on  the  receiver  aperture,  was  43%  [56].  However,  this  kind 
of  receiver  is  not  suitable  for  the  reactions  in  which  liquid  is 
present,  or  in  which  the  products  is  need  to  be  isolated  from  the 
oxygen  in  the  atmosphere  [56]. 

5.  Research  on  SPSRs 

Initially,  the  SPSR  has  been  studied  in  the  early  1 980s  in  order  to 
create  a  direct  absorption  central  receiver  which  is  able  to  interface 
with  high  temperature  (>900  °C)  electric  power  and  chemical- 
production  cycles  [57]. 

Martin  and  Vitko  studied  the  falling  solid  particles,  as  a 
pioneering  work,  in  a  solar  central  receiver  at  SNL  [57].  The  heat 
absorption  characteristics  of  particles  have  been  studied  using  a 
radiant  heat  source  [52].  Possible  particle  materials  have  been 
evaluated  with  respect  to  optical  properties  and  structural  stability 
[58,59].  Some  experiments  have  been  made  on  the  dynamics  of  the 
particles  and  on  the  heat  exchange  [60-62].  The  early  work  on  the 
research  of  SPSR  has  been  reviewed  and  summarized  by  Falcone 
[62-64].  More  recent  SPSR  work  included  on-sun  testing  and 
optical  characterization  of  a  2  MWt  SPSR  prototype  capable  of 
achieving  a  temperature  higher  than  300  °C  at  SNL  [65]. 

Computational  models  have  been  built  to  simulate  receiver 
operation  and  to  aid  in  design  efforts  [51,64].  Many  theoretical 
simulations  have  been  performed  on  the  internal  fluid  dynamics  of 
the  receivers,  in  various  configurations  [53,66-69].  All  these 
studies  simulate  with  global  models  the  radiation  absorption  of  the 
particle  curtain,  the  thermal  transport  between  the  particles,  the 
air  and  the  cavity,  and  the  fluid  dynamics  of  the  air  in  the  cavity; 
these  are  all  fundamental  aspects  to  assess  the  efficiency  and  the 
working  temperature  of  the  plant.  The  particle  temperature  and 
the  efficiency  of  energy  absorption  of  the  receiver  have  been 


investigated  by  Evans  [66-68]  utilizing  numerical  models.  Using 
computational  fluid  dynamics  (CFD)  technique,  the  performance  of 
a  solar  falling  particle  reactor  in  a  pilot  plant  was  investigated  by 
Meier  [53].  The  models,  however,  were  limited  to  two-dimensional 
cases  due  to  a  limitation  in  computational  resources  and  most  of 
these  models  did  not  consider  the  interaction  between  the 
particles  and  the  flow  inside  the  SPSR.  The  3D  numerical  model 
was  built  by  Chen  [69]  and  the  impact  of  the  particle  size  and  other 
operating  conditions  on  the  performance  of  the  SPSR  have  been 
studied  numerically.  However,  the  3D  study  was  confined  in  the 
receiver  only.  In  fact,  the  heat  change  between  the  air  inside  the 
cavity  and  the  air  outside  the  cavity  is  rather  important,  especially 
when  a  wind  field  is  present. 

Windows  for  solar  receivers  are  usually  expensive,  brittle, 
require  careful  mounting,  but  often  fail  to  withstand  the  high  solar 
fluxes  and  high  temperature,  reducing  the  solar  energy  absorption 
efficiency  [56].  Aerowindows  have  the  potential  of  increasing  the 
efficiency  of  an  SPSR  by  reducing  convective  losses  from  an  open 
receiver  aperture  and  eliminating  reflection,  convection  and 
reradiation  losses  from  a  comparable  glass  window.  Aerodynamic 
windows  consist  of  a  transparent  gas  stream,  which  is  injected 
from  an  air  jet  across  the  receiver  aperture  to  isolate  its  interior 
from  the  surrounding  atmosphere  [2,70-74].  Even  though,  the 
wind  conditions  may  still  have  an  important  effect  on  the 
performance  of  SPSRs.  The  aerodynamics  and  thermal  interactions 
between  gas  flow  and  solid  particles  are  very  complicated  due  to 
the  combination  of  convective  and  radiative  exchange.  The 
complexity  increases  considering  the  combined  influence  of  the 
aerowindow  flow  and  the  wind  [2,72-74].  The  wind  effect  on  the 
performance  of  an  SPSR  has  been  investigated  recently,  with  or 
without  the  protection  of  an  aerowindow  [2,72-75].  In  this  present 
design,  the  SPSR  is  only  used  as  a  heat  absorbing  device. 

6.  Solid  particle  selection 

The  technical  feasibility  study  of  solid  particles  has  been 
conducted  and  concluded  by  Hurby  in  SNL  [51,52,61,62,64].  The 
candidate  particle  materials  for  a  SPSR  must  have  at  least  three 
important  characteristics:  (1)  low  propensity  towards  agglomera¬ 
tion,  (2)  high  fracture  resistance,  and  (3)  large  solar  absorptivities. 
Furthermore,  the  particles  must  be  inexpensive  and  readily 
available  in  large  quantities. 

The  initial  material  evaluation  studies  considered  some 
traditional  high  temperature  ceramic  materials  available  in  the 
particle  sizes  of  interest  in  SNL  [51 ,52,61 ,62,64].  Techniques  have 
been  established  for  evaluation  the  agglomeration,  fracture  and 
optical  properties  characteristics  of  the  SPSR  candidate  solid 
materials,  from  alumina,  silica,  silicon  carbide,  zircon,  etc.  Custom 
experimental  devices  for  high  temperature  radiative  optical 
properties  measurements  have  been  designed.  The  agglomeration 
investigation  has  indicated  that  the  best  behaving  pure  ceramics 
for  use  at  1000-1200  °C  are  alumina,  silica  and  zircon.  Fracture 
resistance  studies  have  indicated  that  critical  impact  velocity  can 
be  related  to  the  compressive  strength  and  facture  toughness  of 
the  particle.  It  was  found  that  impact-induced  fracture  is  probably 
not  a  problem  for  small  rounded  ceramic  particles  in  an  SPSR.  The 
effect  of  solar  cycling  on  the  compressive  strength  can  be 
examined  by  cycling  paced  particle  samples  in  a  solar  furnace. 
Toughness  can  be  examined  as  a  function  of  static  heat 
treatment.  Bulk  solar  absorptivity  measure  has  indicated  that 
the  candidate  material  which  performed  the  best  with  respect  to 
agglomeration  had  low  solar  absorptivities  [51 ,52,61 ,62,64].  This 
indicated  that  additional  materials  must  be  considered  in  order  to 
find  a  better  solution  of  the  particle  selection,  which  satisfies  the 
physical  and  economic  characteristics,  in  the  further  develop¬ 
ment  of  SPSRs. 
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7.  Physical  model  for  the  numerical  simulation  of  performance 
evaluation  of  an  SPSR 

A  conceptual  design  of  a  solid  particle  receiver  has  been 
proposed  for  experiments  at  SNL,  as  a  part  of  the  STCH  project 
sponsored  by  the  DOE  in  the  U.S.  [2,69,72-74].  The  working  solid 
particle  curtain  used  in  SNL  for  testing  is  called  CARBO-HSP,  a 
commercial  product  from  CARBO  Ceramics,  which  contains  mainly 
bauxite  and  about  7%  iron  oxide.  It  is  assumed  that  the  particles  are 
spherical  in  shape  (sphericity  of  0.9).  The  detailed  parametric  study 
has  been  conducted  on  this  conceptual  design  and  the  results  will 
be  present  in  the  following  sections. 

The  typical  aerodynamic  and  thermal  processes  in  a  SPSR  are 
analyzed  and  shown  in  Fig.  2.  Solar  energy  is  concentrated  and 
directed  to  the  receiver  from  the  aperture  by  heliostats  on  the 
ground.  The  solid  particles  are  fed  from  the  top  of  the  receiver.  In 
order  to  separate  the  hot  hair  within  the  receiver  and  the  cold  air 
out  of  the  receiver,  an  air  jet  may  be  added  on  the  top  of  the 
aperture.  The  aerodynamic  behaviors  of  the  particles  are  impacted 
by  the  wind,  the  air  jet  flow,  the  air  flow  inside  the  receiver  induced 
by  the  falling  of  particles,  and  the  temperature  difference.  Also,  the 
thermal  interactions  which  include  particle-particle  radiation, 
particle-wall  radiation,  particle-air  convection,  and  air-wall 
convection,  are  needed  to  be  investigated  [2,72-74]. 


governing  equations  for  conservation  of  gas  phase  mass,  momen¬ 
tum,  and  energy  will  be  solved  [2,69,75]. 

A  realizable  k-s  model  is  used  for  closing  the  turbulent  Navier- 
Stokes  equation  system.  A  bibliographical  study  concerning  the 
validation  of  the  CFD  model,  including  the  turbulence  model  (the 
realizable  k-s  model),  showed  that  this  approach  has  been 
extensively  validated  for  a  wide  range  of  flows,  including  the 
channel  and  layer  flows  [2,60,69,75].  Compared  with  the  standard 
k-s  model,  the  realizable  k-s  model  is  likely  to  provide  superior 
performance  for  flows  involving  rotation,  boundary  layers  under 
strong  adverse  pressure  gradients,  separation  and  recirculation 
(including  flows  over  obstacles).  For  all  of  these  cases,  the 
performance  of  this  model  has  been  found  to  be  substantially 
better  than  that  of  the  other  turbulence  models.  The  realizable  k-s 
model  was  also  proved  to  be  a  preferable  model  for  successful  3D 
simulation  of  air  flow  inside  the  SPSR  and  interactions  of  solid 
particles  and  the  air  by  Chen  [69].  The  details  of  the  realizable  k-s 
model  have  been  given  Refs.  [2,69,75],  together  with  the 
differences  with  the  standard  k-s  model  and  the  advantages  of 
the  realizable  k-s  model. 

By  integrating  the  force  balance  on  the  particle,  the  trajectory  of 
the  particles  can  be  predicted  in  a  Lagragian  reference  frame.  This 
force  balance  equates  the  particle  inertia  with  the  forces  acting  on 
the  particle,  and  can  be  written  as 


8.  Mathematical  model 

The  air  flow  inside  the  solid  particle  receiver  can  be  influenced 
by  drag  forces  and  buoyant  convective  heat  transfer  in  the  vicinity 
of  the  particles  because  of  the  interaction  with  the  particles.  The 
particle  volume  fraction  in  the  SPSR  has  been  found  to  be  less  than 
0.1%  [60,69].  Therefore,  the  collisions  between  particles  are 
negligible.  Hence,  the  two-way  coupled  Euler-Lagrange  method, 
which  includes  the  exchange  of  heat  and  momentum  between  the 
gas  phase  and  the  solid  phase,  is  used  to  simulate  the  gas-particle 
flow  in  this  paper.  Particle-particle  collisions  and  the  resulting 
momentum  transfer  are  not  included  at  present. 

With  the  typical  operating  conditions  of  the  receiver,  the 
characteristic  Reynolds  numbers  within  and  out  of  the  receiver  are 
calculated  to  be  larger  than  1.0x1 05.  Therefore,  a  turbulent  flow  is 
expected  within  and  outside  of  the  receiver.  The  time-averaged 


Fig.  2.  Typical  aerodynamic  and  thermal  processes  in  a  solid  particle  receiver  with 
air  jet  flow. 
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where  the  relative  Reynolds  number  Rep  can  be  calculated  by, 
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The  suitable  empirical  equation  for  the  drag  coefficient  for  the 
particles  in  the  solar  receiver  is  proposed  by  Clift  [76], 


CD=^-(-i+0A5Re2p'3)  (3) 

The  effect  of  the  pressure  gradient  and  particle  acceleration  on 
drag  coefficient  can  be  neglected  since  the  material  density  ratio 
between  particle  and  air  is  larger  than  103.  The  energy  equation  for 
the  particle  is  given  by 


PpCp*d3p  ^  =  NuXjtdp(T  -  T„ )  +  n -d2£pcr(T4R  -  Tp)  (4) 

In  order  to  predict  the  interaction  of  an  incident  radiation  field 
with  the  falling  particle  curtain  and  receiver  wall  accurately,  the 
discrete  ordinate  method  with  a  solar  ray  tracing  model  is  used  to 
calculate  the  radiation  field  inside  the  receiver.  The  radiative 
transfer  governing  equation  for  the  spectral  intensity  4(r,  s)  of  an 
absorbing,  emitting,  and  scattering  medium,  at  position  r  in  the 
direction  s  can  be  written  as  [69,75], 

V  ■  (lx(?,s)s)  +  (ax  +  crs)Ix(r,s)  =  axn2lbx 

-y  p4jt 

+  ^yo  (5) 

The  radiative  transfer  equation  is  solved  for  a  finite  number  of 
discrete  solid  angles,  each  associated  with  the  vector  direction  s 
fixed  in  the  Cartesian  system  (X,  Y,  Z)  in  the  discrete  ordinates 
radiation  model.  The  discrete  model  transforms  the  radiative 
equation  into  a  transport  equation  for  radiation  intensity  in  spatial 
coordinates.  The  solar  ray  tracing  algorithm  is  employed  to  predict 
the  direct  solar  illumination  energy  source  that  results  from 
incident  solar  radiation  [69,75].  This  algorithm  takes  a  beam  that  is 
modeled  using  the  sun  position  vector  and  illumination  para¬ 
meters,  applies  it  to  any  or  all  wall  or  inlet/outlet  boundary  zones 
specified,  performs  a  face-by-face  shading  analysis  to  determine 
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well-defined  shadows  on  all  boundary  faces  and  interior  walls,  and 
computes  the  heat  flux  on  the  boundary  faces  that  results  from  the 
incident  radiation  [75]. 

9.  Numerical  procedure 

The  above  governing  equations  are  solved  in  the  Cartesian 
coordinate  system  with  a  control-volume  finite  difference  method. 
A  non-staggered  grid  storage  scheme  is  adapted  to  define  the 
discrete  control  volumes.  A  power  law  differencing  scheme  is 
employed  in  the  present  paper  and  the  solver  used  is  a  segregated 
solver.  In  order  to  resolve  the  coupling  between  pressure  and 
velocity,  the  SIMPLE  algorithm  is  adopted.  The  governing 
equations,  which  were  discrete  and  non-linear,  are  linearized 
using  an  implicit  technique  with  respect  to  set  of  dependent 
variables.  The  algebraic  equations  are  solved  iteratively  using  an 
additive  correction  multi  grid  method  with  a  Gauss-Seidel 
relaxation  procedure. 

9.2.  Wind  speed  classification  and  wind  attacking  direction 

The  wind  speed  classification  based  on  Beaufort-values  is  used 
in  the  present  investigation  [77,78].  The  wind  speed  range  from 
“light  air”  to  “fresh  breeze”  (Beaufort  number  1-5)  will  be  studied. 
With  each  Beaufort  number,  a  typical  speed  is  employed  for  the 
present  model.  The  typical  values  and  the  wind  speed  classification 
are  listed  and  analyzed  in  Ref.  [2].  The  statistics  of  wind  speed  from 
a  local  station  in  Las  Vegas,  Nevada,  U.S.  [79]  are  illustrated  in 
Fig.  3.  From  Fig.  3,  the  selected  wind  speeds  for  studying  (from 
Class  1  to  6)  are  ensured  to  cover  most  of  the  wind  speeds  in  the 
local  region  [2]. 

Fig.  4  shows  an  SPSR  in  a  wind  field  with  an  attacking  angle.  Two 
kinds  of  wind  fields  will  be  investigated.  The  first  kind  of  wind  filed 
is  wind  blowing  parallel  to  the  ground.  The  attacking  angle  a  (in 
plane  of  XOZ),  from  the  top  view  of  the  domain,  stands  for  the  angle 
between  the  wind  velocity  and  the  normal  vector  of  the  aperture. 
The  typical  wind  directions  studied  in  the  present  model  are  a  =  0°, 
a  =  45°,  and  a  =  90°.  Assuming  the  aperture  is  facing  the  south 
direction,  the  wind  is  attacking  from  the  southwest  direction  in  the 
cases  of  a  =  45°.  The  results  with  wind  attacking  from  the 
southeast  direction  {a  =  135°)  are  assumed  to  be  the  symmetrical 
to  the  results  of  a  =  45°.  In  the  case  with  a  =  0°  the  wind  attacking 


Fig.  3.  The  statistics  of  wind  speed  in  Las  Vegas  and  the  typical  values  for  study  in 
the  present  model. 
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Fig.  4.  The  domain  for  CFD  simulation  of  the  wind  effect  on  the  solid  particle 
receiver. 

direction  is  normal  to  the  aperture.  While,  wind  blows  tangential 
to  the  aperture  in  the  cases  of  a  =  90°.  The  attacking  angle  /3  (in 
plane  of  YOZ  in  Fig.  4)  is  used  to  distinguish  the  different  angle  of 
wind  from  the  side  view  of  the  domain.  An  attacking  angle  of 
/3  =  90°  stands  for  wind  blowing  from  the  top  and  /3  =  -90°  stands 
for  wind  blowing  from  the  bottom,  both  tangentially  to  the 
aperture.  An  attacking  angle  of  /3  =  45°  stands  for  wind  blowing 
from  the  top  but  with  an  angel  of  45°  and  /3  =  -45°  stands  for  wind 
blowing  from  the  bottom  and  with  an  angel  of  45°  to  the  vector  of 
the  aperture. 

9.2.  Selection  of  solar  incident  radiation  flux 

The  statistics  of  solar  incident  radiation  from  a  local  station  in 
Las  Vegas,  Nevada,  U.S.  [78].  The  eight  hours  average  values  (from 
8:00  to  16:00)  and  the  ten  hours  average  values  (from  7:00  to 
17:00)  are  calculated  and  plotted  in  Fig.  5.  Typically,  a  uniform 
solar  incident  radiation  flux  of  920  kW/m2  is  employed  for  the 
solar  receiver  performance  calculations  [66,68,69].  In  order  to 
provide  this  amount  of  solar  radiation  flux,  1195  solar  heliostats 


Fig.  5.  The  average  solar  incident  radiation  in  Las  Vegas. 
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units  are  needed  for  the  maximum  values  of  solar  radiation  (in  June 
and  July),  whereas  1559  solar  heliostats  units  are  needed  for  the 
minimum  values  of  solar  radiation  (in  January)  based  on  the  eight 
hours  average  value.  If  ten  hours  operating  time  is  required,  1260 
units  and  1769  units  are  needed  for  the  maximum  and  minimum 
values  of  solar  radiation  in  order  to  provide  solar  incident  radiation 
flux  of  920,920  l<W/m2  [2]. 

9.3.  Geometry 

Different  designs  of  the  SPSR  have  been  proposed  in  SNL.  For 
example,  a  commercial  design  with  5  m  wide,  4  m  deep  and  6.7  m 
high  has  been  proposed  [70].  However,  the  experimental  and 
numerical  data  for  this  commercial  size  receiver  is  still  sparse. 
Another  simple  designer  for  the  distribution  of  the  particles  after 
dropping  has  been  set  up  in  SNL  [70].  In  this  simple  design,  there  is 
no  top  board  and  front  board  and  the  size  is  0.4  m  x  0.4  m  x3  m. 
The  solar  radiation  energy  is  not  load  on  the  particles  and  the 
thermal  aspects  have  not  been  taken  into  accounted.  Another 
conceptual  design  of  the  SPSR  in  SNL  for  laboratory  testing  and 
numerical  investigation  is  shown  in  Fig.  6,  including  the  schematic 
geometry  and  dimensions.  The  geometry  of  this  conceptual  design 
consists  of  a  three  dimensional  rectangular  cavity  that  is  2  m  wide, 
1.58  m  deep  and  3  m  high.  The  width  of  square  aperture  which 
allows  for  the  incident  solar  flux  is  1.5  m.  An  opening  of  0.25  m  by 
1 .5  m  at  the  top  is  used  for  the  introduction  of  particles.  The  heated 
solid  particles  will  be  storage  in  the  catch  hopper  for  the  high 
temperature  solids  storage  after  falling  through  the  receiver  from 
the  top  inlet.  The  parametric  study  will  be  conducted  with  this 
conceptual  design. 

9.4.  Calculating  domain  selection 

The  research  of  wind  effects  on  an  object  is  an  important  subject 
because  of  its  technological  significance,  such  as  for  a  skyscraper  or 


a  bridge  in  a  strong  wind  field.  However,  most  of  the  researchers 
selected  a  much  larger  domain  than  the  object  for  studying.  The 
dimensions  of  the  calculated  domains  are  chosen  randomly,  as 
long  they  are  assumed  large  enough  [80-83].  For  the  present  3D 
numerical  model  considering  the  wind  effect,  the  mesh  for 
calculation  should  be  refined  in  the  receiver,  and  should  be  even 
finer  in  the  region  close  to  the  particle  curtain.  The  mesh  for  the 
wind  field  can  be  coarser.  In  order  to  make  the  mesh  system 
smoother  and  to  save  the  calculating  time,  the  calculating  domain 
is  better  to  be  as  small  as  possible  as  long  it  can  capture  the  main 
phenomena  in  the  receiver.  Therefore,  a  domain  independent 
study  is  necessary.  To  simplify  the  domain  independent  study,  the 
space  dimensions  for  each  side  (up,  down,  left  and  right)  will 
employ  a  same  value,  S,  while  the  approaching  length  is  assumed 
to  be  L,  as  shown  in  Fig.  2.  During  the  domain  independent  study, 
the  approaching  length  3  m,  5  m,  7  m,  10  m,  and  15  m  are  studied. 
Later  on,  the  side  space  is  studied  for  S  =  0  m  (without  side  space), 
1  m,  3  m,  5  m,  7  m,  and  10  m.  As  reported  in  Ref.  [2],  the 
approaching  length  L  =  7  m  and  side  space  S  =  3  m  can  be 
considered  a  good  choice  for  calculation  since  both  the  pressure 
and  velocity  difference  are  less  than  5%  with  more  increase  of  the 
size  of  the  calculating  domain  [2]. 

9.5.  Mesh  independent  study 

After  determination  of  the  calculating  domain,  the  grid 
independence  is  checked.  Four  mesh  systems  have  been  created 
for  numerical  simulation.  Mesh  #1  has  222,939  nodes  and  209,320 
hexahedral  cells;  mesh  #2  has  597,350  nodes  and  624,120 
hexahedral  cells;  mesh  #3  has  1,068,926  nodes  and  1,031,076 
hexahedral  cells;  mesh  #4  has  1,563,626  nodes  and  1,501,236 
hexahedral  cells.  The  differences  for  both  values  decrease  to  less 
than  5%  when  the  mesh  nodes  increases  to  1 ,068,926  (mesh  #3).  To 
save  the  calculating  time  and  to  ensure  the  accuracy,  a  mesh  with 
1,068,926  nodes  and  1,031,076  hexahedral  cells  (mesh  3)  is 
selected  for  the  present  work  [2,72]. 

10.  Numerical  results  and  discussions 

A  two  band  isotropic  radiation  calculation  is  made  with  a  solar 
band  from  0.2  to  1.5  jutm  and  an  infrared  band  from  1.5  to  30  p,m. 
etailed  information  related  to  the  input  conditions  of  the  analysis 
and  the  physical  properties  of  the  solid  particle  are  given  in  Table  1 . 
Diffusely  scattering  bounding  walls  with  an  emissivity  of  0.8  are 
used.  It  is  assumed  that  the  walls  are  perfectly  insulated  to  be 
adiabatic. 

10.1.  The  particle  size  influence 

A  smaller  size  of  the  particles  leads  to  a  higher  cavity  with  the 
same  operating  conditions.  This  conclusion  has  been  demonstrated 
in  Ref.  [69]  with  consideration  with  or  without  the  outlet  opening 
at  the  bottom  of  the  receiver.  Chen  [69]  studied  different  size  of 
solid  particles  from  200  to  1000  jam  and  the  cavity  efficiency  keeps 


Table  1 

Simulation  conditions  and  solid  particle  properties. 


Gas  flow 

Solid  particle 

Operating  pressure:  101, 000  Pa 

Air  outlet  condition:  pressure  outlet 
Wind  temperature:  300  K 

Air  jet  temperature:  300  K 

Back  flow  temperature:  300  K 

Particle  density:  3200  kg/m3 

Heat  capacity:  1085J/kgK 

Thermal  conductivity:  6.67  W/mK 
Particle  emissivity:  0.9 

Particle  diameter:  650  |xm 

Particle  total  flow  rate:  5.0  kg/s 
Particle  initial  temperature:  873  K 
Particle  initial  velocity:  0.088  m/s 
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dropping.  The  same  conclusion  has  been  drawn  by  Chen  [84]  in  her 
simulation  with  consideration  of  cooling  effect  of  the  atmosphere. 
However,  the  solid  particle  size  selection  should  take  into  account 
the  manufacturing  issues  and  a  too  small  size  of  the  particle  may 
cause  the  escape  of  the  particles  from  the  aperture.  Obviously,  the 
particle  distribution  at  the  exiting  opening  is  impacted  by  the 
particle  size. 

2  0.2.  The  influence  of  the  bottom  opening 

The  heated  solid  particle  will  be  released  form  the  bottom  of  the 
receiver.  One  scheme  is  to  keep  the  exit  open  all  the  time  and 
another  scheme  is  to  keep  the  exit  close  and  only  open  it  after 
collecting  a  stack  of  particles.  In  the  case  of  the  receiver  with  a 
bottom  opening,  the  hot  air  exits  with  the  particle  flow  through  the 
bottom  of  the  receiver  when  the  outlet  opening  is  open.  In  the  case 
where  particles  are  allowed  to  exit  the  bottom  of  the  receiver  the 
average  receiver  temperature  is  lower  than  the  averaged  receiver 
temperature  in  the  case  where  the  particles  are  not  allowed  to  exit, 
due  to  the  increased  backflow  of  air  from  the  aperture  [69]. 
Therefore,  a  case  with  bottom  opening  has  lower  cavity  efficiency 
than  a  case  with  a  closed  bottom  opening. 

10.3.  The  wind  effect  without  protection  of  an  aerowindow 

The  cases  with  different  wind  speed  and  attacking  direction  has 
been  examined  without  the  aerowindow  protection.  At  first,  the 
cases  with  different  wind  speed  and  attacking  normal  to  the 
aperture  (a  =  0°,  f$  =  0°)  is  examined  without  the  aerowindow 
protection.  Then,  the  wind  attacking  direction  effect  on  the 
performance  of  a  solid  particle  solar  receiver  has  been  studied  with 
different  angle  of  a  and  different  wind  speeds.  When  the  angle  a  is 
studied,  the  angle  f$  is  kept  to  be  zero  in  all  the  cases. 

Fig.  7  shows  the  cavity  efficiency  for  each  case,  and  Fig.  8  shows 
the  particle  exiting  temperature  for  all  the  cases.  The  results  are 
compared  with  results  from  the  cases  with  normal  wind  velocity 
(a  =  0°).  From  Fig.  8,  the  solid  particles  can  be  heated  up  to  1050  I< 
or  a  higher  temperature  successfully  with  this  small  conceptual 
experimental  design.  An  SPSR  in  a  field  without  wind  undergoes 
large  convection  loss  and  has  very  low  cavity  efficiency  without 
protection  of  an  aerowindow.  When  wind  attacks  normally  to  the 
aperture,  a  stronger  wind  blowing  normal  to  the  aperture  will  keep 


Fig.  7.  The  cavity  efficiency  for  different  wind  speed  and  attacking  angle  a  (/3  =  0°). 


Fig.  8.  The  cavity  efficiency  for  different  wind  speed  and  attacking  angle  ct(fi  =  0°). 

more  hot  air  in  the  cavity  so  that  higher  cavity  efficiency  can  be 
obtained.  The  cavity  efficiency  of  the  cases  with  wind  an  attacking 
angle  a  =  45°  has  no  significant  difference  from  the  cases  with 
normal  wind  attacking.  The  wind  with  a  low  velocity  can  reduce 
the  large  convective  heat  loss  so  the  SPSR  has  slightly  higher  cavity 
efficiency  than  the  normal  wind  cases.  With  a  high  wind  speed 
(larger  than  4.5  m/s),  the  wind  with  an  attacking  angle  a  =  45° 
causes  strong  turbulence  in  the  cavity  and  results  in  a  lower  cavity 
efficiency  than  the  cases  with  a  =  0°.  The  efficiency  and  the  particle 
exiting  temperature  of  the  cases  with  wind  an  attacking  angle 
a  =  90°  keep  the  same  tendency  with  the  cases  with  the  cases  of 
a  =  45°  and  a  =  0°.  Basically,  the  cavity  efficiency  tends  to  increase 
when  wind  spend  increases. 

Then,  the  wind  effect  of  different  attacking  angle  f$  is  studied, 
while  the  angle  a  is  kept  to  be  zero.  The  cavity  efficiency  of  cases  in 
which  wind  blows  from  the  top  (with  positive  /3,  typically  f$  =  90° 
and  f$  =  45°)  is  shown  in  Fig.  9,  compared  with  the  cases  in  which 
the  wind  blows  normal  to  the  aperture.  The  efficiency  of  an  SPSR  in 
a  wind  filed  with  f$  =  45°  increases  greatly  when  wind  velocity 


Fig.  9.  The  cavity  efficiency  for  different  wind  speed  and  attacking  direction  (with 
a  =  0°  and  positive  ft). 


T.  Tan,  Y.  Chen / Renewable  and  Sustainable  Energy  Reviews  14  (2010)  265-276 


273 


increases  from  Vwind  =  2.5  m/s  to  Vwind  =  4.5  m/s  and  then  drops 
slowly  because  the  blowing  wind  can  not  prevent  the  large  loss  of 
the  heated  air  from  the  aperture  at  a  low  wind  velocity.  However, 
the  blowing  wind  serves  as  an  new  “aerowindow”  and  keeps  most 
of  the  hot  air  in  the  cavity,  when  the  wind  velocity  reaches  4.5  m/s 
or  higher.  On  the  other  hand,  a  higher  wind  velocity  will  bring 
more  cold  air  into  the  cavity  and  the  cavity  efficiency  starts  to  drop 
again  after  the  wind  formed  “aerowindow”.  The  cavity  efficiency  of 
the  cases  with  /3  =  90°  is  relatively  low  and  has  no  obvious 
difference  with  the  cases  with  /3  =  0°. 

The  cavity  efficiency  of  cases  in  which  wind  blows  from  the 
bottom  (with  positive  /3,  typically  /3  =  -90°  and  /3  =  45°)  is  shown  in 
Fig.  10,  compared  with  the  case  in  which  the  wind  blows  normal  to 
the  aperture.  From  the  comparison,  the  wind  effect  on  the  cavity 
efficiency  is  more  complex  when  wind  blows  from  the  bottom.  For 
an  SPSR  in  wind  filed  with  /3  =  90°,  the  cavity  efficiency  increases 
with  wind  speed,  reaches  a  maximum  value  at  Vwind  =  4.5  m/s  and 
Vwind  =  6-5  m/s,  and  then  decreases  again.  The  cavity  efficiency 
increases  with  wind  speed  gradually  for  the  cases  with  /3  =  45° 
because  the  wind  blow  from  /3  =  -45°  can  not  formed  an 
“aerowindow”  until  the  wind  speed  reaches  9  m/s.  A  higher  wind 
speed,  however,  reduces  the  loss  of  hot  air  from  the  aperture  and 
let  less  cold  air  to  enter  the  cavity.  As  a  result,  the  cavity  efficiency 
increases  with  the  wind  speed  till  the  “aerowindow”  is  formed  at 
^wind  —  9  m/s  [2]. 

2  0.4.  The  wind  effect  with  protection  of  an  aerowindow 

In  order  to  check  the  protection  of  the  aerowindow  formed  by 
the  air  jet  flow,  cases  with  air  jet  velocity  of  3  m/s,  5  m/s  and  8  m/s 
are  studied.  The  results  are  compared  with  the  case  without  the 
aerowindow  (l/aj  =  0m/s).  Fig.  11  shows  the  cavity  efficiency  for 
each  case.  From  the  comparison,  the  cases  with  air  jet  velocity  of 
8  m/s  have  relatively  higher  cavity  efficiency  than  the  cases  with 
lower  air  jet  velocities,  since  a  low  air  jet  velocity  may  not  able  to 
form  a  protective  curtain  (aerowindow)  to  prevent  the  cold  air 
from  entering  the  receiver.  Obviously,  the  cold  air  blown  in  will 
decrease  the  particle  exiting  temperature  and  the  cavity  efficiency 
by  convection  heat  transfer  [2]. 

For  cases  of  no  air  jet  flow(VflJ-  =  0  m/s),  a  stronger  wind  blowing 
normal  to  the  aperture  will  keep  more  hot  air  in  the  cavity  so  that  a 
higher  cavity  efficiency  can  be  obtained.  The  case  with  no 
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Fig.  10.  The  cavity  efficiency  for  different  wind  speed  and  attacking  direction  (with 
a  =  0°  and  negative  /3). 


Fig.  11.  The  cavity  efficiency  for  different  wind  speed  and  air  jet  velocity. 


aerowindow  and  no  wind  undergoes  a  large  convection  loss, 
which  causes  very  low  cavity  efficiency.  For  cases  of  air  jet  velocity 
Vaj  =  3  m/s,  the  heated  air  in  the  cavity  escapes  quickly  when  there 
is  no  wind.  When  the  wind  speed  increases  to  class  one,  such  as 
1  m/s,  the  blowing  wind  is  so  weak  that  it  is  unable  to  prevent  the 
great  convection  loss  by  the  escape  of  the  heated  air  from  the 
cavity.  Moreover,  the  wind  brings  more  cold  air  into  the  cavity. 
Therefore,  the  cavity  efficiency  decreases.  When  the  wind  speed  is 
increased  to  2.5  m/s,  the  stronger  wind  can  prevent  most  of  the 
escape  of  the  heated  air.  As  a  result,  the  cavity  efficiency  increases 
from  50%  to  54%.  As  can  be  anticipated,  a  higher  wind  speed  will 
more  efficiently  keep  the  hot  air  inside  the  cavity.  After  the  cavity 
efficiency  reaches  a  maximum  value,  it  will  decrease  again.  This  is 
because  a  much  stronger  wind  will  bring  cold  air  into  the  cavity 
and  decrease  the  efficiency  after  most  of  the  hot  air  is  kept  in  the 
cavity.  For  a  relatively  higher  air  jet  velocity,  such  as  5  m/s  or  8  m/s, 
heat  loss  by  escape  of  heated  air  is  high  only  when  there  is  no  wind. 
So,  the  cavity  efficiency  is  increased  when  the  wind  velocity  is 
increased  to  1  m/s.  However,  when  the  wind  speed  keeps 
increasing,  the  cavity  efficiency  decrease  again  since  too  strong 
a  wind  speed  will  cause  cold  air  to  enter  the  cavity.  Therefore,  the 
wind  effect  on  the  performance  of  the  solid  particle  receiver  is  very 
complex,  especially  when  considering  different  air  jet  velocities. 
The  protection  of  the  air  curtain  formed  by  the  air  jet  can  enhance 
the  cavity  efficiency  and  it  is  better  to  keep  the  air  jet  injection 
velocity  not  smaller  than  5  m/s.  From  Fig.  1 1 ,  the  efficiency  is  in  the 
range  from  46%  to  57.  Therefore,  the  cases  of  air  jet  velocity  of  8  m/ 
s  will  be  studied  as  typical  cases,  with  aerowindow  comparing  with 
the  case  without  the  aerowindow  [2]. 

The  wind  attacking  direction  effect  on  the  performance  of  a 
solid  particle  solar  receiver  has  been  studied  with  different  wind 
speeds  and  directions.  The  air  jet  velocity  is  set  at  8  m/s.  The  angle 
a  (see  Fig.  4)  is  first  studied,  while  the  angle  /3  is  kept  to  be  zero. 
Fig.  12  shows  the  cavity  efficiency,  and  Fig.  13  shows  the  average 
particle  exiting  temperature  for  all  the  cases.  The  results  are 
compared  with  results  from  the  cases  with  normal  wind  velocity 
(o'  =  0°)  and  results  without  protection  of  an  aerowindow.  The 
combination  effect  of  the  wind  and  the  air  jet  flow  is  complex. 
Basically,  the  protection  of  the  aerowindow  ensures  the  receiver  to 
have  higher  efficiency.  However,  the  combination  of  the  high  speed 
wind  and  the  air  jet  flow  may  lead  to  high  turbulence  of  the  air 
inside  the  cavity  and  cause  low  cavity  efficiency,  such  as  the  cases 
with  Vwind  >  4.5  m/s  and  a  =  45°.  In  an  SPSR  in  a  wind  field  with 
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Fig.  12.  The  cavity  efficiency  for  different  wind  speed  and  attacking  direction. 


Fig.  13.  The  average  exiting  temperature  for  different  wind  speed  and  attacking 
direction. 

a  =  45°,  the  wind  forces  a  great  amount  of  cold  air  to  enter  the 
cavity  from  the  left  side  of  the  aperture  and  cause  a  great  amount  of 
hot  air  to  escape  the  cavity  from  the  right  side  of  the  aperture.  The 
cavity  efficiency  of  the  cases  with  wind  attacking  angle  a  =  90° 
tends  to  decrease  as  wind  speed  increases,  with  the  protection  of 
an  aerowindow  formed  by  an  air  jet  (Vaj  =  8  m/s).  This  tendency 
coincides  with  that  of  the  cases  of  a  =  45°  and  a  =  0°  but  the 
influence  of  wind  on  the  particles  is  less  and  a  higher  cavity 
efficiency  can  be  obtained  than  the  cases  of  a  =  45°. 

The  wind  effect  of  different  attacking  angle  /3  (see  Fig.  4)  is  then 
studied,  while  the  angle  a  is  kept  to  be  zero.  The  cavity  efficiency  of 
cases  in  which  wind  blows  from  the  top  (with  positive  /3,  typically 
/3  =  90°  and  /3  =  45°)  is  shown  in  Fig.  14,  compared  with  the  cases  in 
which  the  wind  blows  normal  to  the  aperture.  The  protection  of  an 
aerowindow  is  very  helpful  for  the  wind  blow  from  top  (with 
positive  /3).  A  cavity  of  efficiency  54%  or  higher  can  be  obtained  for 
all  the  cases.  For  cases  with  =  45°  and  /3  =  90°,  the  air  jet  flow 
forms  an  ideal  aerowindow  easily  with  the  help  of  the  blowing 


Fig.  14.  The  cavity  efficiency  for  different  wind  speed  and  attacking  direction  (with 
a  =  0°  and  positive  /3). 


Fig.  15.  The  cavity  efficiency  for  different  wind  speed  and  attacking  direction  (with 
a  =  0°  and  negative  /3). 

wind.  A  slight  decreasing  of  the  cavity  efficiency  with  the 
increasing  wind  speed  is  cased  by  the  entrance  of  more  cold  air 
jet  flow. 

Fig.  1 5  shows  the  cavity  efficiency  of  cases  in  which  wind  blows 
from  the  bottom  (with  negative  /3,  typically  =  -90°  and 
/3  =  -45°).  The  wind  effect  is  more  complex  if  the  wind  blows 
from  the  bottom  (with  negative  /3).  High  cavity  efficiency  can  be 
obtained  with  the  protection  of  aerowindow  a  low  wind  speed  for 
P  =  -90°.  When  the  wind  speed  reaches  2.5  m/s,  the  strong  wind 
breaks  the  aerowindow  from  the  bottom  of  the  aperture  and 
causes  cold  air  to  enter  the  cavity,  and  leads  to  a  minimum  cavity 
efficiency.  When  the  wind  speed  reaches  4.5  m/s,  the  strong  wind 
blowing  upwards  generates  a  circle  of  air  at  the  aperture.  The 
cavity  efficiency  is  increased  because  a  part  of  cold  air  enters  the 
cavity  is  sucked  out  again.  This  circulation  becomes  strongest 
when  Vwind  =  9  m/s  and  a  new  aerowindow  can  be  formed.  The 
cavity  efficiency,  consequently,  keeps  increasing  with  the  increas¬ 
ing  of  the  wind  speed  from  4.5  m/s  to  9  m/s. 
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The  wind  with  speed  of  Vwind  =  1  m/s  and  /3  =  -45°  presses  on 
the  air  jet  flow  and  forms  an  ideal  aerowindow  and  high  cavity 
efficiency  can  be  obtained.  When  the  wind  speed  reaches  2.5  m/s,  a 
small  amount  hot  air  starts  to  escape  from  the  cavity  and  the  cavity 
efficiency  starts  to  decrease.  The  cavity  efficiency  reaches  a 
minimum  value  when  the  wind  speed  reaches  4.5  m/s  because 
most  of  the  air  jet  flow  is  blew  away.  A  large  amount  of  hot  air 
escapes  from  the  cavity  and  a  large  amount  of  cold  air  enters  the 
cavity.  A  wind  blowing  from  /3  =  -45°  can  form  an  aerowindow 
from  down  side  when  the  wind  speed  reaches  6.5  m/s  and  9  m/s. 
The  cavity  efficiency  is  enhanced  significantly  with  this  reversed 
aerowindow  [2]. 

20.5.  The  influence  of  the  air  jet  temperature 

The  temperature  of  the  air  jet  flow  can  be  increased  by  a  bypass 
heat  exchanger  if  necessary.  In  order  to  check  how  this  parameter 
affects  the  cavity  efficiency,  two  cases  with  different  air  jet 
temperature  are  compared.  Both  cases  are  under  the  conditions  of 
Vaj  =  8  m/s,  Vwind  =  4.5  m/s.  The  basic  case  has  an  air  jet  tempera¬ 
ture  of  Taj  =  300  K,  and  the  air  jet  temperature  for  another  case  is 
enhanced  to  Taj  =  873  K.  The  particle  exiting  temperature  is 
increased  from  1084.13  K  to  1087.73  K  and  the  cavity  efficiency 
is  enhanced  from  55.82%  to  56.79%  when  the  air  jet  temperature  is 
increased  from  300  I<  to  873  K.  From  the  comparison,  the  increase 
of  the  air  jet  temperature  does  not  enhance  the  cavity  efficiency  too 
much. 

It  was  reported  that  the  cavity  of  the  case  with  Taj  =  873  K  has  a 
larger  medium  temperature  region  compared  with  the  case  with 
Taj  =  300  I<  [72].  This  region  with  medium  temperature  will  cause 
the  cold  air  entering  from  the  aperture  be  heated  up  earlier  and  one 
more  air  circulation  occurs  right  near  the  exiting  region  of  the  solid 
particles.  Besides,  the  air  jet  flow  is  under  a  high  buoyancy  and 
trying  to  flow  up  since  the  air  jet  flow  has  a  relatively  high 
temperature  in  the  case  with  Taj  =  873  K.  he  buoyancy  force 
hinders  the  formation  of  the  aerowindow  and  more  cold  air  enters 
the  aperture,  which  reduces  the  cavity  efficiency  and  decrease  the 
exiting  temperature  of  the  solid  particles  [72]. 

11.  Conclusions 

The  study  on  SPSRs  has  been  reviewed,  including  the  idea 
originality,  the  experimental  work  and  numerical  studies,  and  the 
milestones  achievements.  The  detailed  studies  of  a  US  DOE  funded 
project  on  the  SPSR  have  been  reviewed.  The  review  included  the 
objectives  of  the  STCH  project,  the  conceptual  design  of  an  SPSR, 
the  geometry,  the  particle  properties,  the  calculating  domain 
selection,  the  wind  effect,  the  protection  of  an  aerowindow  and 
other  factors  which  influence  the  cavity  efficiency 

A  smaller  size  of  the  particles  leads  to  a  higher  cavity  with  the 
same  operating  conditions.  However,  the  solid  particle  size 
selection  should  take  into  account  the  manufacturing  issues  and 
a  too  small  size  of  the  particle  may  cause  the  escape  of  the  particles 
from  the  aperture.  The  statistics  of  the  wind  speed  and  the  solar 
incident  radiation  flux  at  a  local  position  are  analyzed.  The 
independence  of  the  calculating  domain  is  studied  in  order  to 
select  a  proper  domain  for  the  numerical  simulation.  In  order  to 
investigate  the  influences  of  the  wind  speed  and  wind  blowing 
direction  on  the  performance  of  the  receiver,  different  wind 
conditions  of  and  different  air  jet  injection  conditions  are 
simulated  numerically.  The  cavity  thermal  efficiencies  are 
calculated  and  the  optimal  injection  conditions  are  analyzed  for 
different  wind  conditions.  From  the  simulation,  the  influence  of  the 
wind  on  the  performance  of  the  SPSR  is  very  complex.  Basically,  to 
form  an  aerowindow  in  front  of  the  aperture  of  the  SPSR  is  helpful 
in  enhancing  the  cavity  efficiency  and  the  particle  exiting 


temperature.  However,  the  combination  effect  of  the  wind  and 
the  air  jet  flow  may  cause  the  efficiency  to  decline  if  the  wind  speed 
is  larger  than  1  m/s.  The  increase  of  the  air  jet  temperature  does  not 
enhance  the  cavity  efficiency  too  much. 
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